Journal of Systematics and Evolution

Understanding evolutionary relationships in hyperdiverse plant groups remains a major challenge in systematics. The orchid genus Bulbophyllum, the second largest genus of flowering plants, represents an exceptional example of phylogenetic and morphological complexity. Relationships, particularly within the species-rich Asian clade, have remained poorly resolved due to extensive morphological variation and limited resolution in previous phylogenetic studies. Here, we reconstructed phylogenetic relationships using 63 plastid genes from 355 specimens representing 322 species and 65 of the 97 recognised sections of Bulbophyllum. Our analyses confirmed that the genus comprises five major evolutionary lineages comprised of species predominantly from Australasia, Madagascar, Continental Africa, Neotropics, and Asia. We provide the first robust phylogenetic evidence for a dichotomous split within the Asian clade into two well-supported lineages: the Asian-Malesian clade and the Malesian-Papuasian clade, with the latter containing a strongly supported Papuasian subclade. Additionally, this study supports the monophyly of several currently recognised sections while clarifying relationships in previously problematic groups. This study provides the most comprehensive plastid-based phylogenomic framework for Bulbophyllum to date and establishes a foundation for future taxonomic revision and integrative analyses of diversification and trait evolution within this hyperdiverse genus.

Although genetic diversity is a fundamental component of biodiversity, we lack data for a majority of species, particularly in biodiversity hotspots such as tropical forests. We present a comparative genetic dataset of 19 tropical tree species (including one palm) from the Caribbean island of Puerto Rico and neighboring islands (Hispanola and the US Virgin Islands). Using a reduced-representation sequencing technique (SLAF-seq), we identified species-specific single-nucleotide polymorphism (SNP) datasets with 24,413 to 433,637 high-quality SNPs per species. The focal species represent a range of life-history and climate associations, which may be relevant to their genetic structure. Therefore, we also include complementary information on species functional traits (wood density, leaf thickness, specific leaf area, maximum height, and seed dry mass), as well as geographic distributions and climatic associations from species distribution models.

The remarkable diversity of life on Earth results from evolutionary processes functioning across different spatial and temporal scales. Species diversification occurs through various mechanisms and at widely varying rates, but identifying the conditions that trigger bursts of diversification over short timescales remains a central challenge in evolutionary biology. This difficulty is more pronounced when incomplete lineage sorting (ILS), hybridization, and ongoing gene flow obscure evolutionary relationships and complicate species delimitation. In this study, we investigated the evolutionary history and species boundaries within a group of recently diverged Petunia lineages shaped by pervasive gene flow. We integrated phylogenomic, population genetic, and species delimitation approaches to reconstruct lineage relationships and assess whether these lineages represent distinct species or stages along a speciation continuum. By applying methods that account for both ILS and gene flow, we revealed that most lineages are not fully independent evolutionary units but rather occupy intermediate positions along this continuum. Gene flow played a crucial role during diversification, blurring species boundaries and generating reticulate evolutionary patterns. Our findings demonstrate that traditional phylogenetic trees may oversimplify relationships in such systems, while phylogenetic networks offer a more accurate representation of evolutionary history. Comprehensive and integrative analyses, such as those employed here, are essential for capturing these complex dynamics. Ultimately, only four lineages could be confidently recognized as distinct species, whereas the remaining represent cases of ongoing divergence. These results emphasize the need to refine species delimitation frameworks for systems characterized by recent divergence and extensive reticulation.

Grasslands, as one of the most important terrestrial ecosystems globally, have root functional traits that serve as key indicators of plant responses to environmental changes and hold significant ecological importance. To reveal the current status, research hotspots, and frontier trends in the field of grassland root functional traits, this study analyzed relevant literature from the Web of Science Core Collection database between 2000 and 2025. It employs bibliometric methods and utilizes visualization tools such as CiteSpace for a systematic analysis. The results indicate that research in this field has been continuously increasing since 2000, reflecting a growing research interest. China, the United States, and Germany are the leading countries in terms of publication output. However, collaboration networks among authors, institutions, and countries are still not tight enough to form a truly global cooperative network. Co-occurrence analysis of keywords and literature clustering reveal that the research hotspots in this field are mainly concentrated in six directions: multidimensional characteristics of root functional traits, interactions between root functional traits and climate change, synergistic effects of root functional traits and soil microorganisms, responses of root functional traits to land-use changes, coupling of root functional traits with ecosystem functions, and applications of root functional traits in agriculture and ecological restoration. Future research should focus on promoting innovation and standardization of research methods, conducting long-term monitoring, deeply exploring the mechanisms of root-microbe interactions, implementing cross-scale integrative research and model construction, and building international collaborative networks.

Animal naming is fundamental to scientific communication, yet it also reflects the historical and cultural contexts in which names are bestowed. Scientific names function as taxonomic labels and enduring records of human engagement with nature. Owing to this dual role, species names have recently attracted increasing attention from historical and humanities perspectives, both for their informative value and for the biases they may encode. To objectively assess these patterns at a large scale, we investigated etymological trends across Animalia using a comprehensive dataset of species names. Our analyses reveal that naming practices are shaped by a combination of historical events, taxonomic traditions, and cultural influences. Major global disturbances coincided with marked declines in species descriptions, whereas advances in biological techniques were associated with shifts in naming practices. Furthermore, etymological trends differed among phyla, indicating that taxonomic communities vary in their naming conventions. These differences suggest that taxonomists preferences, shared aesthetics, available knowledge, and cultural biases are differentially preserved in scientific names. Together, our results demonstrate that zoological nomenclature constitutes a valuable archive for understanding the historical and cultural dimensions of taxonomy.

We present a genome assembly from a specimen of Quercus canariensis (Fagaceae; Fagales; Magnoliopsida). The assembly was generated using PacBio HiFi long reads with an approximate sequencing depth of 39X and scaffolded using a reference-guided approach. The genome sequence has a total length of 816.0 megabases for haplotype 1 and 804.8 megabases for haplotype 2. The two haplotypes are each resolved into 12 chromosomal pseudomolecules, with only 3.48% and 1.36% of sequences remaining unplaced in haplotypes 1 and 2, respectively. Assembly completeness is supported by BUSCO scores of 98.3% and 98.2% complete genes for haplotypes 1 and 2, respectively. Structural annotation identified 51,882 and 46,482 protein-coding genes in haplotypes 1 and 2, respectively. This genome assembly provides the first chromosome-scale reference genome for Q. canariensis, laying the base for future genomic and evolutionary studies in this understudied species of the hybridizing white oak species complex. TaxonomyLineage cellular organisms; Eukaryota; Viridiplantae; Streptophyta; Embryophyta; Tracheophyta; Spermatophyta; Magnoliopsida; eudicotyledons; Gunneridae; Pentapetalae; rosids; fabids; Fagales; Fagaceae; Quercus EBI:txid568684 Quercus canariensis Willd. 1809 (Willdenow)

Phylogenetic trees represent the evolutionary histories of taxa and support tasks such as clustering and Tree of Life reconstruction. Many established comparison methods, including the Robinson-Foulds (RF) distance, assume identical taxon sets. A methodological gap remains for trees with distinct but overlapping taxa. Existing approaches either prune non-common leaves, which can discard information, or complete both trees such that they share the same taxa. Completion is more comprehensive, but current methods typically ignore branch lengths, which are essential for identifying evolutionary patterns. This paper introduces k-Nearest Common Leaves (k-NCL), an algorithm for completing rooted phylogenetic trees defined on different but overlapping taxa. The method uses branch lengths and topological characteristics and does not rely on a specific distance measure. The k-NCL algorithm is designed to preserve evolutionary relationships in the trees under comparison. The running time is O(n2), where n is the size of the union of the two leaf sets. Additional properties include preservation of original distances and topology, symmetry, and uniqueness of the completion. Implemented in Python, k-NCL is evaluated on biological datasets of amphibians, birds, mammals, and sharks. Experimental results show that RF combined with k-NCL improves phylogenetic tree clustering performance compared to the RF(+) tree completion approach. Availability and implementationAn open-source implementation of k-NCL in Python and the datasets used in this study are available at https://github.com/tahiri-lab/KNCL.

The combination of target capture sequencing (TCS) with low-coverage whole genome sequencing (lcWGS), an approach known as Hyb-Seq, has allowed the integration of natural history collections into the genomics revolution, transforming biodiversity research. To implement Hyb-Seq, a collection of genomic targets, often nuclear orthologs, is needed to design probes for TCS. In flowering plants, the universal Angiosperms353 probe set has been proven resolutive at multiple evolutionary scales, with caveats. Malpighiales is known to be one of the most challenging flowering plant orders to resolve. Within this order, the clusioid clade ([~]2.2K species, 94 genera, five families) is no exception. To resolve phylogenetic relationships in this recalcitrant clade, we design a custom probe set, the Clusioids626 kit, composed of 39,936 120-mer probes targeting 626 nuclear orthologs ([~]6.6M nucleotides). This probe set includes all Angiosperms353 targets and 273 clusioid-specific ones, carefully chosen taking copy-number, length evenness, and phylo-informativeness into account. We test our probe set on 70 accessions representing all families and tribes in the clusioid clade. On average, 50.4% of TCS reads mapped to our targets, recovering a median of [~]600 orthologs. Relationships for all clusioid families are fully resolved for our nuclear targets. Additionally, 105 plastid coding DNA sequences were retrieved from the lcWGS fraction. A strong cyto-nuclear conflict was detected. The Clusioids626 kit performs better than the universal Angiosperms353 enrichment panel alone. Our kit design workflow can be extended into other lineages for which a universal probe set exists but more resolution is needed.

Diversisporales comprises species with worldwide distribution that produce glomoid, otosporoid, or tricisporoid spores. The recent reorganization of the order by Oehl et al. (2016) recognizes two families, Diversisporaceae and Corymbiglomeraceae, comprising one and five genera, respectively. Several Glomeromycotan specimens collected in northern and southeastern Mexico and in French Polynesian atolls were characterized using both morphological and molecular analyses. Phylogenetic inference revealed that they represent new members in the Diversisporales, supporting the reorganization of the genus Redeckera into three independent lineages: Albocarpum gen. nov., with A. arenaceum sp. nov., A. leptohyphum sp. nov., and A. fulvum comb. nov., Pulvinocarpum pulvinatum gen. et comb. nov., and Redeckera, which retains five species, including R. varelae sp. nov. In addition, we described Melanocarpum mexicanum gen. et sp. nov. and Diversispora papillosa sp. nov. A broader phylogeny, based on eDNA sequences and representative of Diversisporales species, including the newly described taxa, further supported the split of Redeckera and suggested three additional clades likely corresponding to a new family and two new genera, awaiting the discovery of representative morphospecies to be formally described. Using eDNA sequences metadata, the occurrences of the newly described taxa were mapped, allowing to recognize distribution patterns, mostly in the pantropical zone, distinguish widespread and rare species, and suggest possible endemisms. Finally, the coexistence of species forming large sporocarps (A. fulvum and A. leptohyphum) alongside species forming spores in loose aggregates (A. arenaceum), prompted us to propose a possible sporulation dimorphism in Albocarpum, an argument previously raised to explain the nested placement of Corymbiglomus and Paracorymbiglomus within the Redeckera clade.

Renowned for its aromatic fruits and economic importance, the genus Vanilla poses longstanding taxonomic and phylogenetic challenges. Despite recent molecular studies, a comprehensive species tree is lacking, and the evolutionary processes and historical patterns shaping the genus remain poorly understood. We present a new, comprehensive phylogenomic framework for Vanilla, based on 349 low-copy nuclear genes and 76 plastid loci from the Angiosperms353 probe set, which we used to infer evolutionary relationships, assess cyto-nuclear and gene-species tree discordance, and thoroughly investigate its historical distribution and diversification. Sampling 43% of the genus, our framework resolves phylogenetic uncertainties, clarifies major clades, confirms prior hypotheses, and reveals novel placements, including V. planifolia and Vanilla subg. Gondwana. Discordances are primarily driven by incomplete lineage sorting, particularly in the vanillin-producing clade, with evidence of both ancient and recent hybridization, including a natural hybrid from the Yucatan Peninsula. Biogeographic analyses indicate a Guiana Shield origin ([~]30 Mya), Amazonia as a major diversification source, the Andes as a permeable barrier, and Central America as the main diversification sink. This study provides a robust evolutionary framework for Vanilla, supporting taxonomic revisions, comparative trait analyses, and a deeper understanding of the processes shaping this economically and biologically important orchid genus.

If youve ever complained about a species name thats a mouthful--say, the soldier fly Parastratiosphecomyia stratiosphecomyioides or the myxobacterium Myxococcus llanfairpwllgwyngyllgogerychwyrndrobwllllantysiliogogogochensis--youre in very good company. But could the readability of binomial scientific names cause more than complaints? Could it influence how much species are studied and talked about? We examined a random sample of 3,019 species names spanning 29 phyla/divisions. We tested whether name length and reading difficulty are associated with species representation in the scientific literature (measured via literature mentions) and their visibility to the public (measured via Wikipedia pageviews). Both species name traits showed significant negative relationships with literature mentions and Wikipedia reads. Increasing name length from 10 to 30 characters is associated with a 66% decrease in expected mentions and a 65% decrease in Wikipedia reads, while shifting from the most to the least readable name in the dataset corresponds to 53% and 76% decreases. These patterns are consistent with something familiar: the fickleness of human attention, responding to features of the world that are far from rational. While creativity in naming is a cherished part of taxonomy, a touch of orthographic restraint may ultimately benefit both science and the species themselves--especially among understudied uncharismatic taxa.

Abstract/SummaryHere we untangle the systematics of the Asiatic white-eye complex (Zosterops spp.) to better understand the early stages of a recent island radiation. We adopt an integrative approach involving allelic data, genome-scale single nucleotide polymorphisms (SNPs), and museum-based morphometrics coupled with a comprehensive sampling to provide the most holistic understanding of the group to date. The island lineages of Asiatic white-eyes across Indonesia, the Philippines, East Asia, the adjacent oceanic islands of the Western Pacific underwent a deep split separating Zosterops everetti and Z. nigrorum in the Phillippines from a very rapid radiation including Z. japonicus, Z. meyeni, and Z. montanus in the Philippines, Japan, and Indonesia. Z. nigrorum catarmanensis on Camiguin South in the Philippines was found to be nested within Z. montanus and a species limit between Z. nigrorum populations on Panay and Luzon was strongly supported. Phylogenetic splits and population structure were detected within the clade containing Z. japonicus, Z. meyeni, and Z. montanus. A well-supported split separates a northern group including Northern Philippines Z. montanus subspecies, Z. meyeni, and Z. japonicus from the southerly Z. montanus taxa. This creates a paraphyletic Z. montanus. However, based on speciation rates within the broader Asiatic white-eye clade this break likely does not yet represent evolutionarily independent species lineages. Morphological evolution is taking place within the Asiatic white-eyes especially within the robust, large-billed subspecies of Z. japonicus on the oceanic islands of Japan and in the newly identified yellow-morph of Z. montanus on Camiguin South.

Characterized by the earliest use of pottery, the Jomon culture was a unique Neolithic culture that spread throughout the Japanese Archipelago. Previous archaeological evidence suggests that Jomon hunter-gatherers colonized the southernmost islands, the Ryukyu Archipelago, by approximately 7,000 years before present (YBP). However, genetic characteristics of the Ryukyu Jomon population and its contribution to the modern population have not been elucidated yet. In this study, we newly sequenced 273 modern and 25 ancient (6,700-900 YBP) whole genomes collected across the Ryukyu Archipelago. Our analysis demonstrated a genetic differentiation between the Hondo (Japanese mainland) and Ryukyu Jomon, dating back to [~]6,900 YBP. After the divergence from the Hondo Jomon, the Ryukyu Jomon experienced severe bottlenecks, with an effective population size of [~]2,000. Admixture between the Ryukyu Jomon and migrants from the historic Hondo population occurred [~]1,000 YBP, which corresponds to the widespread adoption of iron tools and agriculture in the Central Ryukyus. Different demographic histories between modern Hondo and Ryukyu populations resulted in different rates of Jomon ancestry in these populations. By providing a new perspective on the peopling of the Ryukyu Archipelago, this study significantly enhances our understanding of cultural transitions in the region.

G-protein coupled receptors (GPCRs) are responsible for translating environmental signals of various types into cellular signals. Over 40 thousand GPCR orthologs have been discovered in the supergroup Unikonta, and around 800 genes encode for GPCRs in the human genome. In contrast to this astonishing variety, only a handful of GPCR-related genes have been reported in vascular plants, a major group within land plants. In an attempt to advance our understanding of plant GPCRs as well as their role in plant cellular signaling, here we present comprehensive bioinformatic analysis that includes phylogenetic hypotheses, in silico structural analysis, and tissue distribution of transcripts. Altogether, our work strongly suggests that GCR1 is the sole genuine GPCR expressed in Embriophyta. Finally, we briefly discuss the potential role of GCR1 in root hairs, the tubular outgrowths in root epidermal cells that are involved in nutrient absorption, environmental interaction, and root development.

AO_SCPLOWBSTRACTC_SCPLOWThis study reconstructs the phylogeny of an expansive set of Neotropical leaf beetles in the subfamily Chrysomelinae. From 33 species in the genus Platyphora Gistel, and an additional 37 species representing 16 beetle genera, five genes, three nuclear, and two mitochondrial, were sequenced and used to obtain a well-supported molecular phylogeny using both Bayesian and Maximum Likelihood. The subtribes Chrysomelina and Doryphorina (sensu Daccordi 1982) were monophyletic, while the genus Platyphora was polyphyletic. The genus Leptinotarsa Chevrolat is confirmed to be distinct from Stilodes Chevrolat. Host plant family was recorded for both adults and larvae using direct observations where possible. Ancestral host plant use was reconstructed using Bayesian trait analyses. A complicated history of host plant switches among a restricted set of plant families is revealed: In the paraphyletic Platyphora, one clade that includes Proseicela and Leptinotarsa had two switches from Asclepiadiodeae to Solanaceae, one switch to Moraceae, and one switch to Malpighiaceae, another Platyphora clade had switches between Asteraceae and Rauvolfioideae, and from Rauvolfioideae to Asclepiadiodeae, with other members of the same clade feeding on Boraginaceae and Convolvulaceae. All species included in the clade containing Tritaenia and Stilodes fed on Malpighiaceae, and all species included in the Cosmogramma and Calligrapha clade fed on Malvaceae.

BackgroundPhylogenomic studies frequently reveal widespread gene tree discordance, primarily arising from incomplete lineage sorting (ILS) and hybridization and/or introgression. Disentangling these processes is especially challenging in rapidly radiating lineages. The genus Oryza, with its rapid diversification and multiple genome types, exemplifies this pervasive phylogenetic incongruence. We integrated multiple genomic datasets including whole-genome resequencing, transcriptomes, and published genomes from diploid Oryza species. Concatenation and multispecies coalescent analyses recovered a robust, congruent species tree, placing the FF and GG genome groups as a monophyletic basal clade, followed by successive divergence of the EE, CC, BB, and AA lineages, a topology differing from some prior hypotheses. ResultsTo assess the sources of discordance, we employed a suite of complementary phylogenomic methods. Quantifying introgression via Branch Lengths (QuIBL)-based model comparisons suggested that [~]74.17% of gene tree-species tree discordance is better explained by post-speciation introgression, whereas only [~]15.56% is consistent with ILS alone. Phylogenetic networks (PhyloNet) and allele-sharing statistics (D-statistics, f-branch) corroborated these results, indicating widespread historical introgression both within and between genome groups. Furthermore, genome-wide scans using the fdM statistic localized introgressed genomic regions, which showed reduced interspecific divergence and were enriched for genes involved in stress responses and metabolism. ConclusionsTaken together, our results demonstrate that historical introgression, not ILS, is the dominant force shaping phylogenetic discordance in diploid Oryza. The integrative phylogenomic framework implemented here, which quantifies the contributions of introgression versus ILS and maps the genomic footprint of gene flow, provides a replicable strategy for resolving complex evolutionary histories in other rapidly radiating lineages.

With the availability of full genomes of an ever-increasing number of species, many recent phylogenetic analyses have focused on datasets with thousands of loci. In the presence of incomplete lineage sorting (ILS), many gene trees will be discordant with the species tree, which can then be estimated with a supertree method. In a separate stream of research, the reconstruction of phylogenetic networks, which aim to detect and visualize reticulations in addition to the dominant phylogenetic tree, has become common practice. Level-1 networks have the feature that every node can be interpreted as an ancestor of a subset of the taxa of interest and no two distinct cycles share a common node. TriLoNet is a method that constructs a level-1 network on all taxa from 3-taxa networks, so-called trinets. This approach is similar to modern supertree methods, but the trinets are assigned based on a single sequence alignment. Here we present TriMouNet (Trinet Multilocus Network), which uses the tree topology and branch length distribution in gene trees of a multilocus dataset to infer best-fitting trinets, together with scores quantifying their statistical support. These trinets are then puzzled together into a network on all taxa in a TriLoNet fashion. Experiments on simulated and real datasets show that TriMouNet can identify reticulations with low false positive rate, if the gene trees are accurate. On the other hand, TriLoNet applied to the concatenation of all loci, tends to predict wrong reticulations as the consequence of violations of model assumptions.

We present new distributional records of Argiope spiders in India, based on more than 10,000 digital images of the genus from the region curated by iNaturalist (www.inaturalist.org). Notable range expansions to India are documented for three species: A. chloreides Chrysanthus, 1961, A. mangal Koh, 1991, and A. sector (Forssk[a]l, 1776). Second, previously unrecorded field characters, updated distributional data, and a re-examination of published descriptions of type material, support the resurrection of A. undulata Thorell, 1887 as a valid species, long treated as a synonym of A. pulchella Thorell, 1881. Finally, we report the first in situ photographic records of live specimens of the rarely documented A. caesarea Thorell, 1897 and A. macrochoera Thorell, 1891. These varied findings for a small and conspicuous taxon highlight the value of online community-science platforms for documenting the arachnofauna of a biodiverse region, as well as illustrate the need for continued taxonomic review, even within well-known genera.

In their recent study entitled "Ancient DNA reveals the co-existence of domestic horses, donkeys and their hybrids in the prehistorical northwestern China", Li and colleagues (2026) report the genetic identification of three horses, three donkeys and four first-generation hinny hybrids dating to 400-160 BCE from the Mazongshan jade mining site in northwestern China. While a re-analysis of their ancient DNA sequence data confirms the horse and donkey identifications, it indicates that the four putative hinny specimens were, in fact, donkeys. This revision removes the primary evidence originally shown for the presence of hinnies at this site. Therefore, new data from the Mazongshan bone assemblage are required to support the proposed role of hinny hybrids as integral components of trans-regional trade networks during the Late Warring States and Early Han periods.

Polyploidization is widely recognised as a major driver of plant diversification, with many species persisting as mixed-ploidy systems where multiple cytotypes co-exist. Polyploids are disproportionately represented among invasive species, yet their role in facilitating biological invasions remains poorly understood. Lantana camara, one of the worlds most successful invasive plants, exhibits remarkable cytotype diversity, but the distribution and evolutionary relationships of these cytotypes in its native and invasive ranges have remained unclear. Here, we characterise ploidy variation and assess genetic differentiation among cytotypes in invasive L. camara populations across India. Flow cytometry of more than a thousand individuals reveals that tetraploids overwhelmingly dominate the invasive range, accounting for more than 95% of individuals, while triploids and hexaploids occur at much lower frequencies. Using genome-wide ddRAD-derived SNP markers from diploids, triploids, tetraploids, and hexaploids, we find no genetic differentiation among cytotypes. Instead, individuals of different ploidy levels cluster together across multiple genetic clusters, consistent with recurrent and potentially independent origins of polyploids. These patterns further suggest that L. camara polyploids likely arise via autopolyploid formation. Together, our results establish tetraploidy as the predominant cytotype in Indias invasive populations and reveal a lack of cytotype-specific genetic structure. These findings highlight the need to investigate the ecological advantages of tetraploids and the mechanisms that generate cytotype diversity, key steps toward understanding how polyploidy contributes to the invasive success of this globally important species.